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The formation energies and electronic structure of europium doped zinc oxide has been deter¬ 
mined using DFT and many-body GW methods. In the absence of intrisic defects we find that 
the europium-/ states are located in the ZnO band gap with europium possessing a formal charge 
of 2+. On the other hand, the presence of intrinsic defects in ZnO allows intraband f — f tran¬ 
sitions otherwise forbidden in atomic europium. This result coorroborates with recently observed 
photoluminescence in the visible red region [J]. 

PACS numbers: 


I. INTRODUCTION 

Doping has been widely used to tailor the electronic, 
magnetic and optical properties of semiconductors. Wide 
band-gap semiconductors such as ZnO are attractive for 
ultraviolet light-emitting diodes, lasers and high-power 
photonic applications. In ZnO, rar-earth elements can 
be incorporated in the material and the long life times of 
the excited states allow for an easy realization of popu¬ 
lation inversion with promising applications in optoelec¬ 
tronic applications [ll-Q . ZnO has a large band gap of 
3.4 eV and a high thermal conductivity, enabling new 
electroluminescent devices. 

Channeling experiments 0 indicate, that rare-earth 
elements in ZnO are preferentially incorporated at cation 
sites. More recent photoluminescence (PL) and pho¬ 
toluminescence excitation (PLE) invest^ations on ZnO 
nanoparticles corroborate this hnding [a |^. PL inves¬ 
tigations of Eu-doped ZnO nanoneedles showed sharp 
emission lines from Eu^"*", suggesting the emission arises 
from intra-4/ transitions, in addition to the ZnO inter¬ 
band emission 0. The f to f transitions are forbid¬ 
den in the isolated atom. However the spherical sym¬ 
metry is broken when the impurity is incorporated into 
the ZnO matrix. Besides Eu, optical emission from rare- 
earth orbitals has been achieved in ZnO nanowires im¬ 
planted with erbium and ytterbium Q and thin ZnO 
films doped with erbium, samarium and europium [n^ 
El. The above experimental investigations have sug¬ 
gested explanations for the optical activation of rare- 
earths in ZnO El 13 by means of substitutional hydro¬ 
gen incorporation [l5l| or formation of defect complexes 

SEE El- 

Theoretical investigations using density functional the¬ 
ory (DET) in the generalized-gradient approximation 
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(GGA) have performed for Eu doped ZnO El) Eli but 
were unable to ascertain the origin of the experimentally 
observed emission in ZnO. The main challenge here is the 
correct description of both ZnO band edges and defect 
states. It is common understanding that the use of local 
exchange-correlation functionals wrongly described the 
ZnO band gap, which could lead to misleading conclu¬ 
sions on the location of the impurity rare-earth /-states. 
It has been shown that the correct description of the band 
gap is of paramount importance for the understanding of 
impurity and defect states in semiconductors [20l - [^ . Be¬ 
sides, intrinsic defects may also play an important role, 
since during ion implantation they are often introduced 
in ZnO. Jiang et al. have shown that for many open 
d- or /-shell compounds the GW approach can provide a 
consistent and accurate treatment for both localized and 
itinerant states. 

Using DET calculations and the GW technique, we 
have previously shown that a complex containing a sin¬ 
gle oxygen interstitial defect and an europium atom sub¬ 
stituting a zinc atom is a probable candidate to explain 
the observed emission in the red region of the spectrum 
in europium doped ZnO nanowires, El- Note that GW 
calculations do not include excitonic electron-hole inter¬ 
actions. 

In this work we show GWq calculations, as imple¬ 
mented in VASP, for Eu-doped ZnO. Besides the work 
reported in Ref. El we have considered the presence of 
the most common defects in ZnO, such as oxygen and 
zinc vacancies and interstitials. The wave functions are 
kept fixed to the GGA level, whereas the eigenvalues are 
updated in the Green’s function only. Here we show that 
this kind of defect has a relatively low formation energy 
and suggest possible mechanisms for the formation of 
such defect. Besides, we investigate other defects and 
calculate their formation energies and electronic struc¬ 
ture. 
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II. COMPUTATIONAL DETAILS 

In this work we have employed density functional the¬ 
ory (DFT) [H, [23 and many-body GW techniques [ 2 ^ 
to investigate the formation energies and electronic struc¬ 
ture of Eu doped ZnO. The projected augmented wave 
method (PAW) [ 2 ^ has been used as implemented in the 
Vienna Ab-initio Package (VASP) [13, [ 33 . A 3x3x2 
supercell containing 72 atoms with a (2x2x1) k-point 
sampling and a cutoff of 500 eV is used to calculate all 
isolated intrinsic defects and complexes. For EU 2 O 3 we 
have used (2x2x2) k-points. For Eu metal we adopted 
a bcc structure with a (6x6x6) k-point mesh. 

A. Thermodynamic properties 

To verify the thermodynamic stability of the investi¬ 
gated defect complexes, we follow the approach derived 
by van de Walle and Neugebauer[32|. The formation en¬ 
ergy of a neutral defect or impurity is defined as: 

Af = £;fnO-defect _ £;ZnO-bulk _ ^ ^ - Ep), 

i 

( 1 ) 

where Etot(ZnO-defect) is the total energy of a defec¬ 
tive supercell and EtotZnO-bulk is the total energy for 
the supercell of pure ZnO. rii is the number of atoms of 
type i (defects or impurities) that have been added to or 
removed from the supercell and fXi is the corresponding 
chemical potential of each species, q is the charge of the 
defect. The Fermi level Ep is referenced with respect to 
the valence band maximum Ey of the host. 

The lower bound of the chemical potential corresponds 
to the total absence of impurities/defects in the material. 
An upper bound on the chemical potential is given by the 
energy of the elemental bulk phase or other solubility- 
limiting phases. Formation of Eu clusters (metal) can be 
avoided by imposing the europium chemical potential to 
obey to: 

MEu ^ MEu(Eu-bulk) (2) 

where Eu-bulk is the total energy of europium bulk 
in a bcc structure. Formation of europium precipitates 
are avoided by considering the bound for EU 2 O 3 , a very 
stable phase of europium: 

/^Eu ^ /^Eu(Eu 203 ) (^) 

The oxygen and zinc chemical potential are not inde¬ 
pendent, but related by fiznO = MZn+MOi where fizn, ^J-o 
and /iznO are the zinc, oxygen and zinc oxide chemical po¬ 
tentials, respectively. is the formation enthalpy 

of ZnO. Combining Eql2] and [3] and these conditions, we 
obtain the following relations: 


Mo = M02-moi + (4) 

The chemical potential of europium oxide can be 
writen as: 

MEU2O3 = + 2 /iEu-bulk + 3^02-moi ( 5 ) 

where is the enthalpy of formation of eu¬ 

ropium oxide and MEu-buik and M02-moi are the europium 
bulk and oxygen molecule chemical potentials, respec¬ 
tively. This results in 

fiEu < + ^lEu bulk - ^AAFfZ^o (6) 

Here A=0 for oxygen rich and A=1 for oxygen poor 
conditions. 


B. Electronic structure 

The many-body methodology that is underlying the 
GW approximation goe s back to pioneering work by 
Kadanoff and Bavm [3^. l33| and Hedin [2^. In GW- 
based approximations, the screening of the Coulomb in¬ 
teraction, that enters the non-local self-energy operator, 
is calculated microscopically. This description of the 
screened potential IF is a main advantage of many-body 
methods over other approaches, like hybrid functional 
DFT, where screening is only included phenomenologi¬ 
cally. 

The single shot CqWo approximation starting from 
GCA calculations often yields too small band gaps. It 
was suggested [sll, [13 that partially or fully selfconsis- 
tent schemes, in which either Green’s function G {GWq) 
or both the Green’s function and the dielectric matrix 
(GW) are updated can improve the agreement with ex¬ 
periments. We should mention that the implementation 
of selfconsistent GIF or GIFo in VASP has the shortcom¬ 
ing, that only calculations that maintain a quasi-particle 
picture are possible, i.e. satellite peaks can not be ac¬ 
counted for. Even though this means that these are tech¬ 
nically not ’’full” GIFo calculations, we will for simplicity 
adopt the notation to call them GIF and GIFo in the fol¬ 
lowing. For detailed information on the implementation 
we refer the reader to Refs. [sTI. I36l - l40l| ■ 

III. RESULTS 

A. Thermodynamic properties 

It is in principle possible to calculate GW total energies 
by using the Galitskii-Migdal formula (dj- In Ref. [d^ 
it was found that for the electron gas, the total energies 
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can be well described using GW schemes. Alternative 
approaches to tackle the formation energy in connection 
with GW methods have been proposed In VASP, 

the calculation of formation energies at GW level is not 
readily possible. Therefore, we have decided to calculate 
the formation energies at GGA level. This can be justi¬ 
fied taking into account that the GGA calculations are 
used as starting point for more accurate calculations us¬ 
ing the GW method. We would like to mention that we 
have performed PBEO calculations for some defects and 
we note that the location of Eu-/-states are strongly de¬ 
pendent on the amount of HE included in the functional. 

Eor europium oxide, the formation enthalpy is calcu¬ 
lated to be AiZ-f (Eu203 )=-14.43 eV compared with the 
experimental value of -16.51 eV [i^. Eor zinc oxide the 
formation enthalpy is AiZ-f (ZnO)=-2.90 eV, which is in 
good agreement with other GGA values and even 
GGA-bU calculations [1^. This discrepancy is proba¬ 
bly due to the wrong description of the oxygen molecule, 
as has been discussed in Ref.ji^. Using these values, we 
can obtain the binding energy of the defect complexes, 
which are defined as 

V^binding = — 2 _^ Cjf (, G 

i 

where is the formation energy of the defects 

and is the formation energy of isolated defects 

calculated according to Eq. O- The calculated values 
are shown in Table|T] and a negative value means that the 
complex is stable. Galculations for intrinsic defects can 
reproduce well previous calculations reported using local 
functionals [i^, Neutral oxygen vacancies have 

a low formation energy under oxygen poor-conditions. 
On the other hand, zinc vacancies can be formed under 
Zn-poor preparation conditions. Oxygen and zinc inter¬ 
stitials, as expected, have high formation energies, due 
to their size [4^. For Eu doped ZnO at cation site we 
obtain a formation energy of 2.42 eV under 0-rich con¬ 
ditions. Incorporation of Eu at interstitial positions is 
highly unfavorable, due to the strain Eu causes in the 
ZnO, leading to a strongly distorted lattice. Similar re¬ 
sults have been reported in Ref. . Gomplex formation 
with zinc and oxygen vacancies and zinc interstitials have 
a much higher formation energy. This can be understood 
by considering size effects, which causes a large strain in 
the lattice. As the defect complexes are formed in exper¬ 
iments under extreme non-equilibrium conditions (like, 
e.g.,ion-implantation) we find it more appropriate to re¬ 
port the binding energies of the defect complexes and 
hence their stability against dissociation. 

As we can see the most stable neutral defect is under 
0-rich conditions the Eu-Oi(o) complex. It is interesting 
to point out that the formation of neutral oxygen defects 
at interstitial sites in pure ZnO has a high formation en¬ 
ergy under thermodynamic equilibrium. However, as has 
been shown in Ref. (49| the diffusion barrier for this kind 
of defect is relatively low (around 0.2 eV). Therefore, 


once this defect is formed under ion implantation (non¬ 
equilibrium conditions), it can rapidly diffuse in the ma¬ 
terial and form complexes with europium atoms. This 
may explain why the complex Eu-Oi(o) is so stable in 
ZnO. 


TABLE I: Formation energies Ef and binding energies Eb for 
neutral intrinsic defects and defect complexes in ZnO. 


Defect 

Ef (eV) 

Ab (eV) 


O-rich 

O-poor 

O-rich 

O-poor 

Euzn 

2.42 

7.44 



Oi(split) 

2.58 

5.48 



Oi(oct) 

3.38 

6.28 



Vo 

3.20 

0.30 



Vzn 

1.37 

4.27 



Znint 

6.68 

3.78 



Euzn + Oi(split) 

4.52 

14.67 

-0.48 

1.75 

EuZn -1- Oi(oct) 

4.4 

14.55 

-1.40 

0.83 

EuZn -1- Vo 

6.1 

10.45 

0.48 

2.70 

EuZn -1- Vzn 

10.9 

13.80 

7.11 

2.09 

EuZn 2nint 

9.88 

12.78 

0.79 

1.56 


Now let us discuss the charged complexes. In FiglT]we 
show the formation energies of several Eu-complexes as a 
function of the Fermi energy at 0-rich and Zn-rich con¬ 
ditions. The top of the valence band Ev for all defects 
calculations was aligned with the top of the valence band 
of the host ZnO using the averaged electrostatic poten¬ 
tial as described in Ref. . We see that the most stable 
structure at 0-rich conditions is a complex of Eu-Oi(oct) 
in the -I charge state. Although this is the thermody¬ 
namically most stable defect, followed by Eu-subst and 
Eu —Vo(—1), we argue that the fact that the neutral 
Eu-Oi(oct) correponds to the experimental f-f transition 
is very appealing. Obviouly we cannot rule out the co¬ 
existence of these other defects during the preparation 
conditions. Further calculations are needed to clarify 
whether these defects also have transitions in the exper¬ 
imentally observed region. On the other hand at Zn-rich 
the defect Eu-subst is the most stable one for Ep — Ev 
between 0 and 0.6. For values lerger than that the de¬ 
fect Eu-Vo in the -1 charge state is the most stable one. 
It is important to note that under Zn-rich conditions all 
defects have a very high formation energy, which is an 
indication that they are unlikely to form at these condi¬ 
tions. 


B. Electronic properties 

Because GW calculations can be performed using dif¬ 
ferent approximations, the results for the theoretical 
band gap of ZnO has been under debate. Parameters con¬ 
trolling the calculations include the number of bands Hil¬ 
ls^, the exchange-correlation potential for the starting 
wave functionj22j| as well as the use of approximate mod- 
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FIG. 1: (color online) Formation energies of several Eu-complexes as a function of the Fermi energy a) 0-rich and b) Zn-rich 
conditions 


els for the screening, like plasmon pole approximations 

Depending on the starting functional and the details 
of the calculation, values between 2.1 and 3.6 eV are 
obtained for GgWn 1^.13^ li^. Isil - ls^ l . between 2.54- 
3.6 for GWO l39l. I54l456| and between 3.2-4.3 for GW 

mm- 

We start by validating the GWq method for bulk ZnO. 
For this purpose we consider a four-atom wurzite ZnO 
unit cell and employ a 8x8x8 k point sampling with 
an energy cutoff of 500 eV. The resulting band gaps and 
energetic positions of the Zii-3d states (with respect to 
the valence band maximum set at zero) for several levels 
of GW calculations are a) PBE: 0.8 and -5.1 eV, b) PBEO: 

3.2 and-7.3 eV, c) HSE06: 2.5 and-7.1 eV, d) PBE-hGVF: 

4.3 and -7.2eV, e) PBE-hGWo: 3.3 and -7.0 eV. 

For the GWq calculation, a cutoff of 200 eV for the 
response functions y, as well as 1024 bands have been 
employed. It has been shown that a large number of 
bands is necessary to obtain properly converged results 
in earlier one-shot GqWo calculations [5ll - l53| . We specu¬ 


late that the deviations to Ref. are due to the rather 
low number of bands in that publication. The results are 
shown in Fig. O We find a band gap of 3.3 eV in reason¬ 
able agreement with the experimental value of 3.44 eV 
^|, as well as with other all-electron GWo calculations 
53, [55|. It should be pointed out that recently norm 
conserving pseudopotentials were introduced in VASP to 
improve GW calculations [s^. However, as shown in 
Ref. [ 5 ^ , we expect these to only lead to minor modihca- 
tions of the bandstructure of up to 0.2eV, not changing 
the main conclusions of this work. It should be noted 
that that use of GIF-based methods has been able to 
predict the positions of lanthanide /-states in NiO [s^ 
and Ce 02 

The center of the Zn-d orbitals is localted at -7 eV, 
in a gree ment with previous GIFo-calculations for ZnO 
[m, 1^ . [ 53 . We want to point out, that a fully- 
selfconsistent GIF calculation results in a much too high 
quasi-particle gap of 4.2 eV. Again this agrees well with 
recent all-electron calculations that find a quasi-particle 
gap of 4.3 [ 53 . 
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FIG. 2: (color online) Density of states of bulk ZnO calculated 
in the PBE+GWo approximation. The vertical line denotes 
the highest occupied state. 


From here on we show only results for GWq calcula¬ 
tions. It has been shown that the position of defect states 
is not strongly influenced by the number of empty bands 
and the response function cutoff [2^ . The relaxed geom¬ 
etry of Eu at a zinc lattice position without the presence 
of intrinsic defects is shown in Fig. [3l The europium dis¬ 
tance to nearest oxygen atoms are 2.23A and 2.27A for 
in-plane and c-direction, respectively. These values are 
slightly larger than the Zn-0 bond lengths in pure ZnO. 
This means that if Eu is incorporated at a Zn lattice po¬ 
sition, it should not disturb the lattice significantly. We 
would like to point out that in the supercell GWq cal¬ 
culations, the ZnO bandgap is slightly to large, around 
3.6-3.9eV. We attribute this to the reduced fc-point sam¬ 
pling in the supercell calculation. We also several tests 
with scGWo calculations, to also investigate the influence 
of a selfconsistent update of also the wave-functions. In 
these tests, we find that quasi-particle energies are only 
changed by up to 0.1-0.2eV between scGWq and GWq 
calculations, supporting the use of the GWq approxima¬ 
tion in the following. 

In Fig. [4] the GWq electronic density-of-states of sub- 
stitutionally incorporated Eu in this geometry is shown. 
The Eu-/ states lie within the band gap, around 2eV 
above the VBM (valence band maximum), hybridize 
weakly with the Zn-s and 0-p states. As expected for 
substitutional Eu at a Zn lattice position, Eu has a for¬ 
mal charge closer to 2+. The / spin-up orbitals are fully 
occupied, giving a total magnetic moment of 7/iB. Fur- 
themore, the Eu d-states lie deep in the conduction band 
of ZnO, Ri2.5eV above the CBM, leading to an energetic 
difference to the Eu /-states of over 3eV. Therefore, our 
results explain why the 4/^ — 4/®/5d^ optical transition 
at 530 nm (ai 2.3 eV) that is often observed additionally 
to the intra-f transitions for Eu incorporation into other 
host materials can not be found experimentally in Eu 



FIG. 3: (color online) Atomic structure around the europium 
impurity calculated within GGA for substitutional europium 
in ZnO. Blue, red and grey spheres are Eu, O and Zn atoms, 
respectively. 


doped ZnO dim- 



FIG. 4: (color online) Total and projected density of states 
for Eu at Zn site in ZnO calculated within the PBEO-l-GWo 
approximation. The vertical dashed line denotes the highest 
occupied state. Positive (negative) values of the DOS denote 
spin up (down). 

A possible way to modify the oxidation state of Eu 
in the ZnO lattice is therefore a change of its environ¬ 
ment. It is known that durin g io n implantation intrinsic 
defects are likely to form d [l3- therefore investi¬ 
gated Eu doped ZnO in the presence of nearby oxygen 
and zinc vacancies as well as zinc and oxygen intersti¬ 
tials. In the presence of a neutral oxygen vacancy the 
Eu-0 distances for neighboring oxygen atoms remain al¬ 
most unchanged, compared to the case without vacancy, 
and increase slightly to 2.24A. The geometry is shown 
in Fig. [5l The Eu-Zn distance is 3.27Aand relaxation of 
farther neighbors is insignificant. 

The electronic structure for this system is shown in 
Fig. O The states due to the oxygen vacancy lie 1 eV 
above the VBM. This is in good agreement with calcu¬ 
lations for pure ZnO (see e.g. Refs. [11,141,145). The 
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FIG. 5: (color online) Atomic structure around the Eu — Vo 
complex calculated within GGA. Blue, red and gray spheres 
are Eu, O and Zn atoms, respectively. 



FIG. 6: (color online) Total and projected density of states 
for Eu incorporated at Zn position in ZnO next to a oxygen 
vacancy calculated within the PBEO+GWo approximation. 
The vertical line denotes the highest occupied state. Positive 
(negative) values of the DOS denote spin up (down). 


Eu /-states are located now at 2.5 eV above the VBM. 
Moreover, we observe a small splitting of the Eu-/ states 
as well as a small hybridization with the Eu-c? and O- 
p states. The splitting is caused by a shift of one of 
the seven Eu-/ states that we assume to be induced by 
changes in the local environment that affect the interac¬ 
tion matrix elements and therefore the charge distribu¬ 
tion around the impurity. 

Our GWo results are in strong contrast to the find¬ 
ings of Ref. [I^, where the Eu-/ states are energeti¬ 
cally located directly above the Vo states. We attribute 
this discrepancy to the choice of the GGA functional in 
Ref. [l^, predicting a much too narrow band gap. On 
the other hand, the presence of the Vq does not signif¬ 
icantly change the formal charge around the Eu atom. 
Moreover, the positions of the Eu-/ spin down and Eu d- 
states remain practically unchanged. Again this defects 
cannot explain the experimental observed /—/ transition 



FIG. 7: (color online) Atomic structure around the Eu — Zni 
complex calculated within GGA. Blue, red and gray spheres 
are Eu, O and Zn atoms, respectively. 

in ZnO. 

The relaxed geometry of the Eu-|-Zni defect is shown 
in Fig. [71 The bond lengths between the Eu atom and 
the nearest oxygen atoms are 2.19A, 2.19A and 2.24A 
in the in-plane direction, and 2.4lA in the c-direction. 
This is because Eu atom relaxes away from the Zn atom. 
The Eu-Zn distance is 2.72A, while the distances to Zn 
second neighbors are 3.29A. 



FIG. 8: (color online) Density of states of Eu doped ZnO 
in the presence of a nearby zinc interstitial calculated within 
the PBEO-fGWO approximation. The vertical line denotes 
the highest occupied state. Positive (negative) values are the 
spin up (down) components. 

The electronic structure of this system is shown in 
Fig. |8l The location of the Eu /-states lies close to the 
conduction band minimum (CBM), around 2.9eV above 
the VBM. We can infer a weak hybridization with both 
Zn-s and Zn-d states of the Zn interstitial atom. Sim¬ 
ilarly, a very small hybridization with the 0-p states is 
found. Here, the formal charge of the Eu atom remains 
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close to 2+ with the spin-up /-states occupied. 



FIG. 9: (color online) Atomic structure around the 
Eu — complex calculated within GGA. Blue, red and 

gray spheres are Eu, O and Zn atoms, respectively. 



FIG. 10: (color online) Density of states of Eu doped ZnO in 
the presence of a nearby split oxygen interstitial calculated 
within the PBEO-I-GWO approximation. The vertical line de¬ 
notes the highest occupied state. Positive (negative) values 
are the spin up (down) components. 

For the oxygen interstitial plus Eu defect, we have con¬ 
sidered two geometries. The first consists of an Eu atom 
at substitutional Zn site next to an oxygen interstitial in 
dumbbell configuration . This configuration is shown 
in Fig. O In this case the lattice ralaxations are more 
pronounced, The Eu atom is shifted away from the in¬ 
terstitial complex, decreasing the Eu-0 bond lengths to 
2.21 A. The Eu-Zn distance is also slightly increased to 
3.29Aand the distances to the oxygen atoms of the inter¬ 
stitial are 2.3lAand 2.33A, respctively. 

The electronic structure is presented in Fig. [101 We 
observe a splitting of the Eu-/ states that we attribute 
to a combination of the changed symmetry in the local 
environment. This affects both the interaction matrix 
elements as well as their screening, so that different Eu 
states are influences differently. While 6 of the Eu-/ are 
located closer to the VBM, the remaining occupied Eu-/ 


state lies about 1 eV higher in energy. We find contribu¬ 
tions to the localized states from both Oint-p and Ohost 
states. However, also in this case, the Eu atom pos¬ 
sesses a formal charge of 2-\- with all Eu-/ spin-up states 
fully occupied. The oxygen interstitial atom possesses a 
magetic moment of O/is, a further confirmation, that no 
electron is transfered from the Eu atom to the oxygen 
site. 



FIG. 11: (color online) Atomic structure around the Eu — Vzn 
complex within GGA. Blue, red and gray spheres are Eu, O 
and Zn atoms, respectively. 



FIG. 12: (color online) Density of states of Eu doped ZnO 
in the presence of a nearby zinc vacancy calculated within 
the PBEO-I-GWO approximation. The vertical line denotes 
the highest occupied state. Positive (negative) values are the 
spin up (down) components. 


The geometry of Eu doped ZnO in the presence of 
a neutral Zn vacancy (Vzn) is shown in Fig. 1111 The 
Eu atom shifts towards the vacancy, increasing the in¬ 
plane Zn-Eu distance to 3.49A and 3.40A in comparison 
to 3 . 3 A that would correspond to Eu exactly on Zn lat¬ 
tice position. Consequently, the in-plane Eu-0 distances 
also change, asymmetrically, to 2.09A, 2.19A and 2.24A, 
respectively. Along the c-direction the Eu-0 distance is 



























2.22A. The corresponding electronic structure is shown in 
Fig.[T2l We find no occupied Eu-/ states in the band gap. 
These states lie within the VB, close to the 0-p states 
and hybridize both with with 0-p and Zn-d states. The 
local projected moment on the Vzn site is IpB, which 
is aligned anti-parallel to the magnetic moment on the 
Eu which is 6 ^b- However the fact that the Eu-/ states 
are energetically located inside the VB and the unoccu¬ 
pied Eu-/ state is energetically located at 2 eV above the 
VBM, does not fit with experimentally observed transi¬ 
tion. 

The other defect complex including an oxygen inter¬ 
stitial consists of an oxygen at a octahedral interstitial 
site next to an Eu at Zn lattice position. The corre¬ 
sponding atomic structure is shown in Fig. [TSl We find 
much more pronounced lattice relaxations in this case, 
also involving second nearest neighbors. The Eu-Zn dis¬ 
tance is 3 . 2 A, while the Eu-0 bond lengths are 2.2A to 
the oxygen atoms in the basal plane and 2.25A to the 
oxygen along the c-direction. The distance between the 
Eu atom and the oxygen interstitial is 2.24A. This can be 
explained by looking at the coordination Eu adopts. The 
Eu-0 distance to the nearest neighbors is very similar to 
EU 2 O 3 , which possess a coordination number of 6±1 and 
a radial distance of 2.33 ±0.015 A according to Ref[^. 



FIG. 13: (color online) Atomic structure around the Eu-|-0°/t 
complex within GGA. Blue, red and gray spheres are Eu, O 
and Zn atoms, respectively. 

The corresponding electronic structure is shown in 
Fig. [m We find the occupied Eu /-states are located 
in the band gap, with the Eu-/ orbitals being occupied 
with six electrons. The occupied and unoccupied states 
are located at -0.7 eV and 1.3 eV, respectively. We find a 
small hybridization with O-pi(oct) a-nd O-p(host)- The to¬ 
tal magnetic moment of the complex is 5pB, determined 
by a magnetic moment of 6/iB spin up at the europium 
aligned anti-parallel to the magnetic moment of I/tb at 
the oxygen interstitial atom. Therefore, we suggest that 
the intra-4/ luminescence of the ZnO:Eu samples is most 
likely due to Euzn±0°'^* complexes. We have also ex¬ 
plored the mechanism for optical activation of Eu in ZnO 
reported in Ref. [isj ] , where substitutional hydrogen at Zn 
lattice position next to an europium atom is suggested. 
However, this configuration does not change the formal 



FIG. 14: (color online) Density of states of Eu doped ZnO 
in the presence of a nearby octahedral oxygen interstitial cal¬ 
culated within the PBEO-I-GWO approximation. The vertical 
dashed line denotes the highest occupied state. Positive (neg¬ 
ative) values are the spin up (down) components. 


charge of Eu. 

In order to summarize the results discussed above we 
show the position of the defect levels with respect to the 
band edges of ZnO. The results are shown in Table [III 
We show the relative position of Eu-f states for the com¬ 
plexes with respect to the top of the valence band of ZnO. 
The position of the Eu-f states were taken approximately 
in the middle of the /-bands, since we cannot describe 
correctly with DFT the multiplets of the f-orbitals. 

TABLE II: Relative position of Eu-f states for the investigated 
complexes with respect to the top of the valence band of ZnO. 
The position of the Eu-f states were taken approximately in 
the middle of the f-bands. All values are given in eV. 


Complex 

Evb — Ef 

Euzn 

2.8 

Euzn ± Oi(split) 

1.5/3.2 

EuZn ± Oi(oct) 

2.5 

EuZn ± Vo 

2.5/3.2 

EuZn ± Vzn 

1.0 

EuZn 2llint 

3.8 


IV. CONCLUSIONS 

In conclusion, we have investigated Eu-doped ZnO us¬ 
ing DFT and the many-body GW technique within the 
GWo approximation. We find that the position and for¬ 
mal charge of the Eu-f states is strongly dependent on 
the environment around the Eu atom. We conclude that 
the optical activity of Eu in ZnO is most likely due to 
EuVOf^* defect complexes, possibly either in a neutral 
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or -1 charge state or in the presence of both. Finally, 
we believe our results can open the pathway for a better 
understanding of these complexes in zinc oxide. 
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